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Abstract. We extend the APES (Amplitude and Phase Estimation) method of spectral analysis 
to the case of non-quasi periodic signals like appears in OFDM. 



Jean-Philippe PREAUX00 
Notations 

Domic un nombrc complcxc z G C on note respectivement z* et \z\ le conjugue et le mod- 
ule de z. On note i G C lc nombrc complexe verifiant i 2 = — 1. Pour une matrice (ou un 
vecteur) M a coefficients complexes on note respectivement M T , M* et M c ses matrices trans- 
posee, transconjuguee et conjuguee complexe (M* — (M C ) T ). Si M a P lignes et Q colonnes, 
ses lignes seront numerotees de a P — 1 et ses colonnes de a Q — 1. On considere la norme 
euclidienne ||M|| = E?=o \M jk ?)^ ■ 

1. LE CAS UNIDIMENSIONNEL 

Soient 3 entiers strictement positifs, M, L, N, avec N = M +L—1. On considere {yo, yi, . . . , yjv-i} 
une suite de donnees complexes a analyser, et pour t = 0, 1, . . . , L— 1, on pose y(t) = (yt, ■ . ■ , yt+M~i) T 
S C M . Soit u> G [0, 27r[ la frequence etudiee ; pour tout entier K > on considere ajf(w) = 
(f , e 4W , . . . , e l ^ _1 ^) T G C A '. L'ecriture s'avere plus concise en introduisant les deux matrices 
M x L a coefficients complexes, 

Y = (y ,yi, • • .,yi-i) 

A(w) = a Jl /(w)a2(o;) 

et Ton souhaite exprimer : 

Y = a(u)A(uj) + res 
ou a(oj) G C designe l'amplitude de frequence uj du signal. 

f .f . Le filtre APES. Soit h(w) G C M le vecteur contenant les coefficients du nitre en frequence 
bj, Le filtre APES (cf. pQ) est la solution h(w) au probleme d'optimisation sous contrainte en 
h(w) G C M et a{uj) G C : 



min || h*(w) Y — a(aj)a^ (w) || avech*(w) ajf (a>) = f 

h(w) ,a(cj) 



soit encore, 



L-l 

min V |h*(w)y(t) - a(uj)e lut \ 2 avec h* (w) a M M = 1 

h(uj),a(a») ^ 

Le filtre h(iu) est choisi de sorte que : 

(i) le signal traite h*(w) Y soit le plus proche possible au sens des moindres carres de la sinusoide 
discrete a(w)a^(w) de frequence ui et d'amplitude complexe a(oj), 

(ii) le filtre ne distort pas la sinusoide discrete aif(w) de frequence u>. 

II admet pour solution : 

a A PEs(u) = h* APES (u)g(w) 
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avec : 



g(a;) = ix!y(t)e-- i = iYa£(a;) 
t=o 

QH = T(£y(*)y*(*)) -gM-g» = |yy* - g ( w ) g » 



t=0 

1.2. APES generalise aux signaux a bruit. On generalise cette methode au cas avec bruit. 
Considerons une suite xq, xi, ■ ■ . , xn-i d'elements dans C, representant le bruit. Pour < K < N 
on note x^ = (xq, x%, . . . , xk-i) t ■ On note X re / la matrice M x L dont l'element ligne i colonne 
j est Xi-\-j. On souhaite trouver l'amplitude complexe verifiant : 

Y = a(uS) X re j o A(w) + res 

oil o designe le produit de Hadamard de matrices. On cherche le filtre qui verihe les deux conditions 
suivantes : 

- (i) Le signal filtre h*(w)Y soit le plus proche possible (au sens des moindres carres) de la si- 
nusoide a bruit a(u)(xj,oai(u)), 

- (ii) le filtre ne deforme pas les sinuso'ides a bruit xjh o ajn(w). 

Trouver un filtre h(w) verifiant les conditions (i) et (ii) nous ramene au probleme d'optimisation 
en a(ui) et h(w) sous contrainte suivant : 

min |h*(u;) Y — a(oj) x^ o a^(a;)|| avec h*(u)(xjfoai/(w)) = l 

soit encore : 

L-1 

min V \h*(u)y(t) - a(uS)x(t)e iut \ 2 avec b*(u) (x M ° a M (w)) = 1 (*) 

" v ' t=0 

Proposition 1 (filtre NAPES). Le probleme d'optimisation (*) ci-dessus a pour solution (lorsqu'elle 
existe) : 



^■N APES ( w 



ajvAPiJs(^) = h^r AP£;s (a;)g(w;) 

Q _1 (w) (xju o a M (w)) 



(x M ° a M (w))* Q _1 (w) (xm o a M (w)) 



avec 



SM = Tj HT2 E y(*) = (x L o a L ( W )) c 

II x l|| 2 ^ I|xl|| 2 



avec ||x L || 2 = J2t=o \ x t\ 2 - 



1.3. Demonstration. Nous suivons dans ses grandes lignes l'argument figurant dans [T] concer- 
nant la methode APES en le generalisant au cas avec bruit c'est a dire a la resolution de (*). 
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Prcuve. Pour plus de lisibilite nous notons a = a(uj), a = a(w), h = h(w), g = g(w) et Q = Q(w). 
Arrangeons la quantite I(a, h) a minimiser : 

L-l 

/(a, h) = l h * y(*) - a 



„iwt |2 



t=0 
L-l 



= £ (h*y(t) -aa^e*" 4 ) (y'Wh-a'x'fr"') 
t=o 

L-l L-l 

= h* £y(i)y*(i)h+£H 2 |z(i)| 2 
t=o t=o 

L-l L-l 

- a*h ^ a;* (i) e - iwt y(i) -aj] z(*)e iwt y* (*) h 



t=o t=o 

On pose R = || x ^||2 J2t=a y(*) y*(*) ct S = commc donnc ci-dessus, 

/(a, h) = ||x L || 2 h* Rh + ||x L || 2 (|a| 2 - a*h* g - ag* h) 

= ||x L || 2 (h*Rh+|a-h*g| 2 -|h*g| 2 ) 

= ||x L || 2 (h*(R-gg*)h+| a -h*g| 2 ) 

se minimise pour a = h* g. II reste a determiner le h optimal. On s'est ramene au probleme 
d'optimisation : 

minh*Qh sous la contrainte h* (xm o a) = 1 

h 

avec Q = R — gg*. II s'agit d'un probleme quadratique avec une contrainte egalitaire lineaire, 
ct il est facile dc verifier en appliquant les multiplicateurs de Lagrange qu'il admet pour solution 
(lorsqu'elle existe) : 

h _ Q 1 (xm ° a) 

(x M o a)* Q-^xm o a) 

ce qui donne la solution optimale recherchee. □ 



2. LE CAS BIDIMENSIONNEL 

Soient 6 entiers strictement positifs M, L, N, M' , L' , N' avec N = M + L- lctN' = M' + 
V — 1. Soit Ynm' une matrice N x N' a coefficients dans C representant une serie de donnees 
2-dimcnsionnclles. Pour t = 0, 1,...,L — 1 et t' = 0, 1,...,L' — 1 soit y tjt < le vectorise de la 
sous-matrice M x M' de Yn,n' debutant ligne t, colonne t' : 

( Y N , N ,(t,t') ••• Y N , N ,(t,t' + M' -1) 



MxM' 



y t ,t> = vec 

\Y NiN .(t + M-l,t) ••• Y N , N ,{t + M-l,t' + M 1 -I) 

C'est un vectcur dc C M+M ' . Soit w £ [0, 2tt[ ct pour tout cnticrs P, P' , 

ap,p/(w,w') = ( 1 e iw ' ••• e '(^'-i)"' ) T ® ( 1 e 4 " ••• e l{p -^ u ) T 

ou ® designe le produit de Kronecker ; ap t p>(u),u)') est un vecteur dc C PxP . 
Un filtre bidimcnsionnel est une matrice M x M' a coefficients complexes. 

2.1. Filtre APES bidimensionnel. C'est le filtre h(u, u') qui avec l'amplitude complexe a(cj, to') 
est solution du probleme d'optimisation: 

L-l L'-l 

„, min , „EE vec*(h( W , W '))y t> f -a(o;,o;')e^ + " Y ) | 2 

h(u,,u,'),a(u,,u,') t=Q t/=Q 



vec*(h(w,w'))a.M,M'(w,u') = 1 
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II admet pour solution : 

a APES (u,uj') = vec*(h(o;,u/))g(w,a/) 

Q _1 (w,cj')a M ,M'(w,cj') 



vec(h APES (uj,uj')) 



avec, 



a M,M'( u i u ')Q 1 (lo)sl M ,M'( CJ ^ CJ ') 
^ L-l L'-l 

^'^l^ee 

t=0 t'=0 

q(c, w ') = -^-7 E E y*.*' y*V - w') g(' ■ ' 



L-l L'-l 
t=0 t'=0 



t=0 t'=0 

2.2. Filtre APES bidimentionnel pour signaux a bruit connu. C'est le nitre h(u>,u>') qui 
avec a{uj,uj') est solution du probleme d'optimisation : 

L-l L'-l 

,r T n , « E E I vec*(h)y tit / -a(w,w')a;(t,t / )e <(h '* +w '*' ) | 2 (**) 

v " v ' ; t=0 t'=0 

vec*(h(cj,cj')) (x M ,M' °a M ,M') = 1 

avec x(t, t') G C representant le bruit en t, t' et xjj m' = {x(0, 0), . . . , x(M — 1,0),..., x(0, M' — 
l),...,x(M-l,M'-l)) T . 

Proposition 2. Le probleme d'optimisation (**) precedent se minimise pour 

a NA PEs{u,u') = vec(h NAPES (co,w'))* g(cj,cj) 

Q _1 (w,w') (x M ,M' oajj, M '(w,w')) 



vec(hjvAP_Es(w,w')) 



avec 



(xm,m' °a M ,M')*Q 1 (w,o; / )(xm,m' ° &m,m' ) 

L-l L'-l 

XL,L'" n 



g( w , w ') = ^-^E E ^(t^Je-^'^y^t') 
IIX^'II t=0 1 , =0 

1 L-lL'-l 

Q(a;,a/) = 2 E E ^OftM') - g(w,a/)g>,w') 



ll x L,L'„ t=Q t<=0 

o« ||x ilL HI 2 = E^, 1 E?=o 1 W*,* , )l 2 - 

2.3. Demonstration. La preuve suit toujours le meme argument. 
Preuve. On note a = a(co,co'), h = h(u,u}') et 3lm,m' — 

L-l L'-l 

E E \yec(hry(t,t')-ax(t,t')e^ t+ ^\ 2 
t=o t'=0 
L-l L'-l 



E E (vec(h)'y(M') -K(M')e' MWt,) ) x (y*(t, t 1 ) vec(h) - a* x*(t, f') e -*MWt')^ 

t=0 t'=0 

L-lL'-l L-lL'-l 

vec(h)* £ E y(t,f)y*(t,f)vec(h) + EE l« | 2 |^(*,^)| 2 

t=0 i'=0 t=0 i'=0 

L-lL'-l L-lL'-l 

a*vec(h) ]T E :c*(t,t / )e- < ( w * +w ' t '>y(t,t / )-a E E zCM'K^+^V^*') vec(h) 



t=0 t'=0 t=0 t'=0 
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On pose R = j^j, Ef-Q 1 y(t, f) y*(i, t'), 

et g = g ( w , w ') = EfJo 1 E^o 1 ^(*.*0e- i(wtWt,) y (*,*') 

= ||xl,l'|| 2 (vcc(h)* Rvcc(h) + |a| 2 - a*vec(h)*g - ag* vec(h)) 
= \\*l,l>\\ 2 (vec(h)*Rvcc(h) + \a - vcc(h)* g| 2 - |vcc(h)*g| 2 ) 
= ||x L , L ,|| 2 (vcc(h)* (R-gg*) vcc(h) + |a - vec(h)* g| 2 ) 

se minimise pour a = vec(h)* g. II reste a determiner le h optimal. On s'est ramene au probleme 
d'optimisation : 

min vcc(h)*Q vcc(h) sous la contrainte vec(h)* (xm,m' ° slm,m') = 1 

h 

avec Q = R — g g* . On obtient : 

Q 1 (xm.at ° &M,M') 



vcc(h) 



(xm,M' ° &M,M')*Q 1 ( X M,M' ° a M.M' ) 

ce qui donne la solution optimale recherchee. □ 



3. NAPES APPLIQUE AUX SIGNAUX PARCELLAIRES 

Soit N un entier non nul ; considcrons une suite de donnees yo, yi, ■ ■ ■ , Vn-i dont certaines sont 
manquantes, que l'on partitionne en ses sous-suites de donnees connues dc longueur N 2 k-i (k = 
1,2,.. .) entrecoupees par ses sous-suites de donnees manquantes de longueur N 2 k (k = 1,2,. . .), de 
sorte que les N\ premieres donnees sont accessibles, les iV 2 suivantes sont manquantes, et ainsi de 
suite. Soit P le nombre de segments disjoints ; ainsi EiT=i ^ k = ^ ■ Notons y a et y u les 2 vecteurs 

de dimensions respectives N\ + N 3 + h Np et N 2 + N4 + h Np_i contenant respectivement 

la suite de donnees connues et la suite de donnees inconnues dans yo,yi, ■ ■ ■ ,yN-i- On considere 
aussi la suite des bruits (x n )o< n <N = xq,x\, . . . , xn-i- 

On applique NAPES pour obtenir une estimation de a(co) et h(w) de la facon suivante : choisir 
une longueur de nitre initiale Mq. Si on verifie la condition suivante, qui nous permet de construire 
une matrice R de covariance de Mq lignes constitucc des donnees accessibles : 

53 max(0, N k - M + 1) > M 

fc=l,3,...,P 

appliqucr l.a, sinon appliquer l.b, ci-dessous. 

l.a. Soicnt L k = N k — Mq + 1, et notons J le sous-ensemble de {1, 3, . . . , P} pour lequel L k > 0. On 
note, si L\ > 0, xi le vecteur dont les coordonnees sont donnees par la sous-suite xq, x\, . . . ,Xl-i 
de (x n ) n debutant au rang et de longueur L\ (si L\ > 0), etc. . . , et pour tout P 6 J, xp sa 

sous-suite debutant au rang EfceJ ^ k et ^ e l° n g ueur Lp. 

Appliqucr NAPES aux segments de donnees correspondant a J a l'aide des redefinitions suivantes 

N 1 +N 2 + ---+N k ^ 1 +L k 

Y{uj) = v — h^p E E y(WK iwt , 

l^kejW fell keJ t =N 1 +N 2 + -+N k _ 1 

N 1 +N 2 + ---+N k _ 1 +L k 

R = v. .h^ E E y(*)y*0 



Eke.J H x fell 2 fceJ 



t=N 1 +N 2 + -.-+N k _ 1 



l.b. Appliquer NAPES avec Mq = N/2 a la suite de donnees yo,Vi, ■ ■ ■ ,Vn-i cn posant que y„ 
est une suite de zeros. 



2. On determine ensuite une estimation de y„ basee sur l'estimation initiale des spectres a(wfe) et 
h(wfc) determinee precedemment. Elle est obtenue a l'aide de la minimisation suivante : 
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K-l L-l 

min £ E l h *K)yW - u(cuk)x t e iujkt \ 2 
y " k=0 t=o 

Considcrons la matrice L x N a coefficients dans C : 



Hfe 



/ h*K) 

o h*K) 



V o 



et 



z k = a(uj k ) x* o 



h*(w fe ) / 
1 \ 



i(L-l)w fc 



Lc critere dc minimisation devient 



K-l 

En h * 



yi 



- Zk\ 



\ Vn-i ) 

Considerons les matrices Ak et Bk telles que 

( Vo \ 



2/i 



A k y a + S fe y„ 



et soit : 

dk — Zk - A k y a 
De sorte que le critere a minimiser en y u devienne : 

K-l 

E \\B k yu-d k \\ 2 

k=0 

minimise pour : 

K-l K-l 

y u = (Y,B* k B k )- 1 (J2 B *kdk) 

k=0 k=0 

3. Une fois qu'un estime de y u a ete trouve, on re-estime le spectre en appliquant NAPES cette 
fois-ci a y et y u . Cela consiste a minimiser en h(wfe) et a(wfe) le critere suivant : 

£ J2 |h> fe )y(i) - a^Mi^'T 
fc=0 t=o 

Qui se scinde en fait en if problcmes d'optimisations distincts pour determiner h(cjfc) et a(wfc), 
pour k variant de a K — 1. 

L'algorithme d'optimisation consiste alors en : 
Soit 8 > la marge d'erreur imposee 

Appliquer l'etape initiale 1 pour estimer h(u>fc) et a(wfe) pour fc = 0, 1, . . . , K — 1 
FAIRE : 

Estimer y u a l'aide de l'etape 2. 
Estimer h(cjfe) et a(u>k) l'aide de l'etape 3. 
TANT QUE Les valeurs ont change de plus que 5. 
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II s'agit en fait d'une optimisation cyclique du probleme suivant : 

, r?u, u EEl h *K)yW-«K)^We^T 

y„,{a(«»),h(u;»)} k=Q t=Q 
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